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Abstract 

The tensor susceptibility of QCD vacuum is calculated in the global 
color symmetry model. The input parameters for gluon propagators 
are determined via the simplified equation for calculating the pion 
decay constant. The reason for the great discrepancy between our 
results and those from QCD sum rules and from chiral constituent 
quark model is discussed. 
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The tensor susceptibility of QCD vacuum, like the quark condensate or 
the gluon condensate, reflects the non-perturbative aspects of the QCD vac- 
uum directly. It is argued that, the tensor susceptibility is related to a 
chiral-odd spin-dependent structure function that can be measured in the 
polarized Drell-Yan process The earlier estimations for the value of 

tensor susceptibility were obtained by QCD sum rules techniques (jS]-|H]) or 
from chiral constituent quark model P^. Two decades calculations show that, 
the global color symmetry model (GCM) JU] describes the nonperturbative 
aspects of strong interaction physics and hadronic phenomena at low ener- 
gies quite well ([llj- 29j), we naturally expect that GCM is applicable in the 
estimation of the tensor susceptibility of QCD vacuum. Recent investigation 
shows that, the value of tensor susceptibility calculated from an effective 
quark-quark interaction is much smaller than the others Questions are 
therefore arising. What is the reason for this discrepancy? Is anything wrong 
with GCM? This letter aims to answer these questions. 

The QCD partition function for massless quarks in Euclidean space can 
be written as 

Z = J VqVqVAe'^^^''^^^^ (1) 

with the QCD action 

S% q.A]= j dx[g(x)(^ - ig4)q{x) + \f;,F;J\, (2) 

where = A^^, F^^, = d^A^ - duA"^ + g^'^^'A^^Al. By introducing the 
functional W[J] defined as 

the QCD partition function can be rewritten as 

Z = J VqVqe~ ■f'^''^^''^^'^'-''^e^^'^^'^''^'^\ (4) 

The functional W[J] has the expansion 

W[J] = \j Myj;{x)Df,{x,y)Jl{y) + Wr[JI (5) 

where = D'^^^{x — y) is the gluon 2-point Green's function, and 

VTiifJ] involves the higher order n{> 3)-point Green's functions. The GCM 
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is obtained through the truncation of the functional W[J] in which only 
Dlll{x,y) is retained. This model maintains global color symmetry of QCD, 
but the local color SU(3) gauge invariance is lost. For simplicity we use a 
Feynman-like gauge — y) = 5^y5°'^D{x — y). The important dynamical 

characteristics of local color symmetry are included in D{x). The exact form 
of D{x) is not well known. Instead, we use a phenomenological gluon propa- 
gator, which is required to exhibit the properties of asymptotic freedom and 
infrared slavery. The justification of such truncation relies on the successes 
of various calculations. 

The partition function of this truncation can be given as ^21 1^ 

Zgcm = J VqVqexp J dxq^q - ^- j dxdyj;ix)D;l{x - y)jt{y)^ , 

with the quark color current j^{x) = qix)'j^^q{x), or equivalently, 

Zgcm = / I?gPgPAe-^GCM[5,q,A] 

with the GCM action 

SGCM[q,q,A] = J dx[q{x){^-zg4)q{x) + J dxdy^A;{x)[D;l{x-y)]-'Al{y). 

(8) 

By the standard bosonization procedure, the resulting expression for the 
partition function in terms of the bilocal field integration is 

Zgcm = lvB'exp[-S[B']), (9) 

where the action is given by 

g 1 , /"j J B%x,y)B%y,x) 

S[B] = -TrLn[G ] + J dxdy .^.^^^ _ . (10) 

and the quark inverse Green's function is defined as 

G-\x, y) = -y) + A'^'(x, y) . (11) 

Here the quantity arises from Fierz reordering of the current- current in- 
teraction term in Eq. © 

aX = ii^Y^Ai^Y^ii (12) 
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and is the direct product of Dirac, flavor SU(2) and color matrices: 



A' = ^(/d,^75, ^7/., ^7m75) ® {^If, -^Tp) ® (^Je, -^A^) . (13) 

Here we consider Np = 2 flavors as in Ref. [Oj. 

The vacuum configurations are defined by minimizing the bilocal action: 
0, which gives 



SS[13] 
SB 



Bo 



B^oix) = g^D{x)tT[Go{x)A'^] . 



(14) 



These configurations provide self-energy dressing of the quarks through the 
definition S(p) = A^-Bg(p) = ifi[A{p^) - 1] + B{p^). According to Ref. [Ej, 
the self-energy functions A and B satisfy the Dyson-Schwinger equations, 



Bip') 



8 f d^q 



3 i (27r)4 ^ ^' 



A{c^^q ■ p 



q^A\q^) + B\q^) ' 
B{q^) 
q^A\q^) + B^{q^) 



(15) 
(16) 



The quark Green's function at Bq is given by 



Go{x,y) = Go{x - y) 



dV -i^A{p^) + B{p^) 

(27r)>2^2(p2) +^2(^2)' 



;ip-{x-y) 



(17) 



The hadron properties follow from considering deviations from these vac- 
uum configurations. If we consider only the isoscalar a{x) and isovector vf (x) 
fields, the approximate local-field effective action can be taken as JU] 

5[(T,7r] = -TiLn [^A{x -y) +m6^^\x -y) + V[a,n]B{x -y)} 

+ d''z[a^{z) + TT^z)] J d^wB{w)tT[G{w)] , (18) 

where V[a, n] = a"(2^) + i75'7f(^|^) ■ r and m is the quark bare mass. Ex- 
panding the spectrum of S[a, vf] to second order about its minimum S'[l, 0] 
with a{x) = 1 + 6{x) 

S[l+6{x), n^)]-S[l, 0] = ^j? J [{d,6)' + d'z+^fj J [{d^n)' + m^'n'] d'z+- 

(19) 
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it is found that 



and 



^TT^ Jo 



sds 



-^^^ B^{s)[B^{s) - sA\s)] 







B\s)[sA\s) - B^{s)] 



2sA^(s) + 



sA\s) + B\s) 



3m 



sds 



A^(s)B^(s 



sds 



5fsl 



SA2(S) +52(s) ' 

B%s) 



2 + 



SA2(S) +52(s) 



(20) 

(21) 
(22) 

(23) 



[sA2(s) +52(s)]2 

In Eqs. (1211) and (j2SI), all those terms involving the derivatives of A{s) and 
B{s) with respect to s are neglected. Let A{s) = 1, -B(s) = M, with M the 
mass of the constituent quark, Eq. reduces to 



^7r2 Jo 



sds 



2 + 



(24) 



When s approaches infinity, the integrand of Eq. behaves like 2M'^/s, 
which reproduce the result of Ref. [Oj strictly. While s — > 0, it behaves like 
3s/M2 rather than 2s/M2 of Ref. p. This difference arises from the second 
term in the brackets of Eq. (j2H). Can this difference bring some serious 
problems? Let put this question aside at present. 
The tensor susceptibility x is defined as 



X 



where 



47r2 



lyo) 

Q{qq) 



sds 



(25) 



Bis] 



SA2(S) +52(s) 

is the quark condensate, we need only to calculate nx(0) defined as pUj : 

1 2 



(26) 



^nv(o) 



47r2 Jo 



In this equation, if we let A{s) 
reproduced. 



sds 
hB(s) -- 



Bis) 



SA2(S) +52(s)J ■ ^^^^ 

: M, the result of Ref. jH] is strictly 
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Before the numerical calculation of tensor susceptibility, we make an anal- 
ysis of the nontrivial solutions A{s) and B{s) to the Dyson- Schwinger equa- 
tions (fT3j) and (fT^ . Because the phenomenological gluon propagators exhibit 
the properties of asymptotic freedom and infrared slavery, when s = in- 
creases from to oo, B{s) (> 0) decreases from some finite non-zero value to 
zero, while A{s) decreases from some non-zero value down to 1. Therefore, 
A{s) is universally greater than 1. The replacement of A{s) by 1 in Eq. (j27|) 
will lead an increase in the value of n^(0)/12, or in other words, the value 
calculated from Eq. ^T7\\ is expected to be smaller than that obtained by 
the others. This is verified in the recent calculation [30], and will be further 
checked in the calculation below. 

The tactics in calculation is similar to that of Ref. [HOj. The input param- 
eters are adjusted to reproduce the pion decay constant in the chiral limit 

= 87 MeV via Eq. 

To give a convincible conclusion, we choose three different gluon propaga- 
tors. The ultraviolet behavior of these model gluon propagators are different 
from that in QCD [221 ESI- They are model 1: 

g'D{s)=A7r'd^^, (28) 
s"^ + A 

model 2: 
and model 3: 

9'D{s) = 4n^d + . (30) 

+ A s ln[s/A"^ + e\ 

Here d = 12/(33 - 2Nf) = 12/29 and A = 200 MeV. For model 1, the 
self-energy functions A{s) and B{s) varying with s are showed respectively 
in Figs, m and m with the input parameters A = 0.1 GeV^, A = 1.780 GeV. 
Obviously, A{s) is not less than 1. In Table ^ the values of n^(0)/12 for 
model 1 are displayed, and the corresponding values for quark condensate 
(qq) are also listed. It should be noted that the values for quark condensate 
are roughly around those obtained from QCD sum rules (see Ref. j3I] and 
references therein), due to the term B (s) / {sA'^ (s) + (s)) in Eq. against 
the term [B{s)/{sA'^{s) + ^^(s))]^ in Eq. (EB . 

The results for model 2 and model 3 are given in Table |21 It is shown that 
the values of the quantity n^(0)/12 are still very small as that of Ref. ^30ij . 
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Table 1: The values of 11^(0)/ 12 for model 1 with Eq. used to calculate 
The quark condensate {qq) is also presented. 



A [GeV^] 


A [GeV] 


-{qqY'^ [MeV] 


n;,(0)/12 [GeV] 


10-1 


1.780 


279 


-0.0017 


10-2 


1.350 


244 


-0.0015 


10-^ 


0.955 


210 


-0.0013 


10-6 


0.770 


196 


-0.0013 



Table 2: The numerical results for models 2 and 3 with Eq. used to 
calculate 





model 2 




model 3 


A [GeV"] 


A [GeV] n^(0)/12 [GeV] 


A [GeV^ 


A [GeV] n^(0)/12 [GeV] 


2.000 


2.94 -0.0021 


10-1 


1.71 -0.0016 


0.200 


1.51 -0.0015 


10-^ 


0.95 -0.0012 


0.020 


1.44 -0.0012 


10-^ 


0.71 -0.0012 



To further check if the second term in the brackets of Eq. (jzHj) give some 
serious modifications to our results, we drop this term intentionally and ob- 
tain: 

3 A%s)B%s) 
At:^Jo [sA\s) + B\s)f ^ > 

Let A{s) = 1, B{s) = M, the result of Ref. jU] is strictly reproduced. With 
Eq. ()31|). similar calculations can be performed for the various models as 
previously did. For example, for models 2 and 3, the results are given in 
Table 01 One can see that, these results do not make much difference to the 
previous ones. 

To summarize, we have calculated the QCD vacuum tensor susceptibility 
based on the modified version of calculating pion decay constant. The various 
calculations in this letter and in Ref. [20] show that, it is a fact that the 
value of tensor susceptibility calculated in GCM is very small. In these 
calculations, the basic characteristic {A{s) > 1) of quark propagator Go{p) or 
quark self-energy determined from the Dyson-Schwinger equations keeps 
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Table 3: The numerical results for models 2 and 3 with Eq. (j31|) used to 
calculate 





model 2 




model 3 


A [GeV"; 


A [GeV] n^(0)/12 [GeV" 


A [GeV^" 


A [GeV] n^(0)/12 [GeV 


2.000 


2.955 -0.0023 


10-1 


1.755 -0.0018 


0.200 


1.600 -0.0017 


10-^ 


1.010 -0.0014 


0.020 


1.570 -0.0015 


10-7 


0.765 -0.0014 



unchanged. This is probably the main reason to the small value of tensor 
susceptibility. So, for the calculation of QCD vacuum tensor susceptibility, 
GCM formulated to date deviates seriously from the QCD sum rules and 
constituent quark model. 
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Figure 1: The self-energy function A{s) as a function of s for the gluon 
propagator g^D{s) = with A = 0.1 GeV^ A = 1.780 GeV 
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Figure 2: The self-energy function B{s) as a function of s for the gluon 



propagator g D(s) 



487rl 

29 s^+A 



, with A = 0.1 GeV^ A = 1.780 GeV 
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